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ABSTRACT 

We present observations of a series of 10 outbursts of pulsed hard X-ray flux from the 
transient 10.6 mHz accreting pulsar GS 1843-02, using the Burst and Transient Source 
Experiment on the Compton Gamma Ray Observatory. These outbursts occurred reg- 
ularly every 242 days, coincident with the ephemeris of the periodic transient GRO 
J1849-03 (Zhang et al. 1996)), which has recently been identified with the SAS 3 source 
2S 1845-024 (Soffitta et al. 1998). Our pulsed detection provides the first clear iden- 
tification of GS 1843-02 with 2S 1845-024. We present a pulse timing analysis which 
shows that the 2S 1845-024 outbursts occur near the periastron passage of the neutron 
star's highly eccentric (e = 0.88 ± 0.01) 242.18 ± 0.01 day period binary orbit about a 
high mass (M c > 7M ) companion. The orbit and transient outburst pattern strongly 
suggest the pulsar is in a binary system with a Be star. Our observations show a long- 
term spin-up trend, with most of the spin-up occurring during the outbursts. From the 
measured spin-up rates and inferred luminosities we conclude that an accretion disk is 
present during the outbursts. 
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1. Introduction 

The pulsar GS 1843-02 was discovered with Ginga 
in 1988 (Makino 1988) during a galactic plane scan 
conducted as part of a search for transient pulsars 
(Koyama et al. 1990b). These observations found 
6 new transient X-ray sources near galactic longi- 
tude I = 30°. Their locations and inferred distance 
of 10 kpc (based on hydrogen column density esti- 
mates) placed them within the "5 kpc arm" , an in- 
ner spiral arm of the galaxy identified with radio, in- 
frared and CO molecular line observations (Hayakawa 
et al. 1977). GS 1843-02 was detected with a 
flux of 8 mCrab. Periodic pulses with frequency of 
10.5432±0.0002 mHz were detected, with a pulse pro- 
file having a deep narrow notch seen in all energy 
bands. A cross scan located the source to a 50' x 6' er- 
ror box (90% confidence). The SAS-3 source 2S 1845- 
024 was within the original Ginga error box (Makino 
1988), but is outside of the refined Ginga error box 
(Koyama et al. 1990a). 

Koyama et al. (1990b) concluded the six new 
sources they had discovered in the 5 kpc arm were all 
transient accreting pulsars in Be-star binary systems. 
Accreting pulsars in binary systems with Be (or Oe) 
stars form the largest class of known accreting pulsars. 
With the exception of X-Perseus, all of these systems 
are transients. Be stars are main-sequence stars of 
spectral type B that show Balmer emission lines (see 
Slettebak 1988 for a review). This line emission and a 
strong infrared excess is associated with the presence 
of a cool circumstellar envelope. There is strong ob- 
servational evidence that this envelope has the form 
of a disk, extending along the Be star equatorial plane 
(e.g. Quirrenbach et al. 1997). In Be/neutron star 
binary systems the transient outbursts are thought 
to be fueled by accretion from this circumstellar en- 
velope, frequently occuring near periastron passage. 
A review of Be/neutron star binary systems is given 
by Apparao (1994). 

Zhang et al. (1996) analysed data from the Comp- 
ton Gamma Ray Observatory (CGRO) spacecraft's 
Burst and Transient Source Experiment (BATSE) for 
transient hard X-ray emission for the region of the 
sky near I — 30° ,b = 0° . Using Earth occultation 
analysis techniques they identified six hard X-ray out- 
bursts between April 1991 and February 1995 from a 
source, designated GRO J1849-03, which was local- 
ized to a 30' x 60' error box (90%) that contained 
both 2S 1845-024 and portions of the GS 1843-02 error 



box. The outbursts were regularly spaced by 241±1 
days, lasted about 13 days, and reached a peak flux of 
75 mCrab (20-100 keV). The outburst behavior sug- 
gested a Be/neutron star binary, with a natural can- 
didate being GS 1843-02. Pulsations from this source 
during the outbursts were searched for in the BATSE 
data, but not detected at that time. Observations 
with the BeppoSAX wide field camera (Soffitta et al. 
1998) during a predicted outburst of GRO J1849-03 
in September of 1996 identified this source with 2S 
1845-024. Zhang et al. and Soffita et al. suggested 
the identification the source with the transient pulsar 
GS 1843-02. Pulsations were, however, not detectable 
in the BeppoSAX data due to limited sensitivity, leav- 
ing this conjecture unconfirmed. 

Here we present CGRO/BATSE observations of 
outbursts of pulsed flux from GS 1843-02 which are 
coincident with the outbursts of GRO J1849-03 = 2S 
1845-024 detected in BATSE Earth occultation mea- 
surements. Observations of pulsations not detected 
by Zhang et al. (1996) were made possible by sensi- 
tivity enhancements resulting from improvements in 
BATSE pulsar analysis techniques (see section 2.1). 
Pulse timing analyses of these observations are pre- 
sented which show that the outbursts occur near the 
periastron passage of a wide, highly eccentric orbit 
about a high mass companion. In addition these 
observations show that the pulsar spins up during 
outbursts. We conclude from the measured spin-up 
rates and inferred mass accretion rates that an accre- 
tion disk is present during the outbursts. We suggest 
that the pulsar's orbit is inclined from its companion's 
equator, with an outburst initiated and an accretion 
disk formed during each passage of the pulsar through 
its companion's circumstellar envelope. 

2. Observations and Data Analysis 

BATSE is an all sky monitor designed to study 
gamma-ray bursts and transient source outbursts in 
the hard X-ray and soft gamma-ray bands (Fishman 
et al. 1989). The observations reported here use data 
from the eight Large Area Detectors (LADs), which 
are 2025 cm 2 area by 1.24 cm thick Nal (Tl) scin- 
tillators operated in the 20 keV to 1.8 MeV band. 
These are located at the eight corners of the CGRO 
spacecraft. The LADs are uncollimated, with each 
detector viewing half of the sky. The large fields-of- 
view result in high background rates in the detectors, 
which are generally not easily modeled or subtracted. 
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The flux of a discrete source can be sampled twice 
per spacecraft orbit by fitting the steps in count rates 
that occur when rises or sets over the Earth's horizon 
(Earth occultation analysis, see Harmon et al. 1992). 
For pulsars with periods shorter than ~ 500 s, the 
variable portion of a pulsar's pulse profile can be de- 
tected by using Fourier analysis, epoch- folding, or by 
fitting count rates over short intervals, using a model 
that includes a description of the pulse profile, and 
a polynomial which accounts for the combined back- 
ground and the average source flux (pulsar analysis, 
see Bildsten et al. 1997). 

Our pulse timing analyses use the DISCLA channel 

1 data, which consists of discriminator rates for the 
20-50 keV band, continuously read out from all eight 
detectors with a resolution of 1.024 s. Results from 
Earth occultation analyses use the 16 energy channel 
CONT data, which is continuously read out from all 
eight detectors with a resolution of 2.048 s. 

2.1. Blind Frequency Searches 

We have made blind searches for pulsations consis- 
tent with the 10.54 mHz pulse frequency of GS 1843- 

02 (Koyama ct al. 1990a) in the BATSE DISCLA 
data from 1991 July 16 to 1997 June 20. Zhang ct 
al. (1996) searched the same data during the six out- 
bursts of GRO J 1849-03 they observed, resulting in no 
clear detection of pulses from GS 1843-02. Through 
a combination of longer integration times, optimized 
combination of rates from different detectors, and im- 
proved handling of systematic effects due to interfer- 
ing sources, our current searches are approximately 

3 times more sensitive then those made by Zhang et 
al. (1996). With these improvements, and the ex- 
tended search interval, we clearly detect 10 outbursts 
from GS 1843-02, at times consistent with the out- 
burst ephemcris of GRO J1849-03. 

The data were divided into four day intervals, with 
each interval searched for pulsations in the barycen- 
tric frequency range of 10.47-10.61 mHz. The increase 
of the search interval from the one day used by Zhang 
et al. (1996) to four days provides a factor of two 
improvement of sensitivity for a constant frequency 
source. No additional increases in the search inter- 
val duration were made because of the increased dan- 
ger that frequency changes due to intrinsic spin-up 
(or spin-down) and orbital motion would result in a 
loss of coherence (and hence sensitivity) in the search 
(Chakrabarty et al. 1997). 



For each four day interval DISCLA channel 1 data 
was selected for which the source was visible, the 
high voltage was on, the spacecraft was outside of 
the South Atlantic Anomaly, and no electron precipi- 
tation events or other anomalies had been flagged by 
the BATSE mission operations personnel. The rates 
were combined over detectors using weights optimized 
for a source with an exponential spectrum dN/dE = 
(A/E)exp(-E/kT) with temperature kT = 25 keV 
(see Bildsten et al. (1997) appendix A). This weight- 
ing provides an improvement of sensitivity over the 
cosine of aspect angle weighting used in Zhang et al. 
(1996). These combined rates were grouped into seg- 
ments of ?s300 s duration. This segment duration 
was chosen because it contains several pulse periods, 
and yet is short enough that a quadratic model gen- 
erally provides a good fit to the detector background. 
Segment boundaries were chosen to avoid inclusion 
of occultation steps from the bright sources Cyg X-l, 
GRO J1744-28, and GRO J1655-40. 

A pulse profile estimate was obtained from each 
segment by fitting the combined rates with a model 
composed of a sixth order Fourier expansion in pulse 
phase (representing the pulse profile) combined with 
a quadratic in time (representing the background plus 
mean source rate). This resulted in an estimated 
pulse profile (represented by 6 Fourier coefficients) 
for each data segment. The phase model used in this 
initial fitting had a constant barycentric frequency of 

= 10.54 mHz. Six harmonics were used so that the 
narrow notch seen in the Ginga pulse profiles could 
be reasonably represented. 

The set of (typically several hundred) segment pro- 
files for a given four day interval were then searched 
for a pulsed signal using a generalization of the Z\ 
test of Buccheri et al.(1983). The Z\ test assumes 
Poisson statistics. However, the noise in our data is 
often dominated by the aperiodic activity of Cygnus 
X-l, which is usually in the field of view. Our test 
statistic, which we call the Y n statistic, accounts for 
this frequency dependent excess noise. 

We represent the pulse profile in the m th 300s data 
segment as 

n 

r(t) = 5ft amh cxp(z2tt^o W) (1) 

h=l 

where denotes the real part, n = 6 harmonics, a m h 
is the estimated complex Fourier coefficient for inter- 
val m and harmonic h, and the initial phase model 
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is <po(t) = (t — t)vq, with t the barycenter corrected 
observation time and the epoch r = JD245000.5. 

We fit the segment profiles in a four day interval 
to a mean profile which is represented as 



r{4>) = 5ft ^ [in exp(i27r/i</>) 



(2) 



h=i 



where Hh are the estimated Fourier coefficients of the 
mean pulse profile. The pulse phase within segment 
m is represented as 



where 



0(t) = fo(*) + A0 



A<j> m = (t m - t)Av 



(3) 



(4) 



with i m the midpoint barycentric time of segment to, 
and Av a frequency correction estimated for the four 
day interval. Here we assume the phase change within 
a 300 s segment caused by Av is a small fraction of a 
cycle. 

The \ 2 of the fit has the form 



X 2 = J h{Vh) 



h=l 



with 



M 



\a m h - exp(z27r/iA</>„ 



(5) 



(6) 



m— 1 



rah 



where M is the number of 300s segments in the four 
day interval and a m h is the Poisson error of either the 
real or imaginary part of a mn . 

The Z 2 statistic may be expressed as 



£ 

h=l 



= ]T [Mo) Jr n ] (7) 



h=l 



were the J^ m 



is the minimum of J n which occurs 
at the mean Fourier coefficient n™™, and a^ h is the 
formal (i.e. Poisson statistical) error of /x™ m . To 
account for the non-Poisson errors, we multiply a 2 h 
by the reduced X 2 factor J,™ in /(2M-2). This results 
in the following test statistic: 



h=i 



2M-2 

jmin 



(8) 
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Fig. 1. — Blind search results for pulsations from GS 
1843-02. For each four day search interval the plot 
shows peaks in the Yq statistic with value above 35. 
The symbol area is proportional to Yq — 35. The ar- 
rows indicate the peak times of the outburst of GRO 
J1849-03 reported by Zhang et al. (1996). For refer- 
ence JD 2448500.5 = 1991 September 1.0. 



For each four-day search interval the Y& statistic 
was calculated for a grid of frequency offsets. The 
results are summarized in Fig. 1. This shows for each 
search interval the peaks in Yq with value above 35, 
with the symbol size in proportion to 16 — 35. Monte- 
Carlo calculations show that with no signal present 
we should expect «75 peaks on the plot with Y@ > 
35, and a 50% chance of one peak with Y$ > 45. 
We actually find 55 with Yq > 45 in Fig. 1, with a 
maximum value of Yq = 140. 

A series of ten outbursts is clearly evident, as is 
a long-term spin-up trend in the pulse frequency of 
w 2.7 x 10~ 13 Hz s _1 . The arrows give the peaks of 
the outbursts of GRO J1849-03 detected by Zhang 
et al. (1996) using BATSE occultation techniques. 
These are consistent with the pulsed outbursts we 
have detected. The regular outbursting is typical 
of many of the Be star/pulsar systems observed by 
BATSE (Bildsten et al. 1997). 
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Fig. 2. — Frequency changes during the outbursts. 
Pulse frequencies were determined by searches over 
a narrow frequency band chosen to include the out- 
bursts seen in Fig. 1. The detected frequencies minus 
a linear trend (see text) are plotted against the ob- 
servation time folded at a 242 day period, beginning 
with MJD 48362, which provides an approximate or- 
bital phase. A highly eccentric binary orbit combined 
with accretion induced spin-up caused the observed 
signature. The periastron epochs given in Table 1 are 
at phase 44 days in this plot. 

Fig. 2 examines how the observed pulse frequency 
changes during the outbursts. For each four day in- 
terval we narrowed our frequency search range to a 
14 fjttz band centered on the linear frequency model 
i'modei(i) = vq + (t — t)z>, with vq = 10.585 mHz, 
t = JD 2450000.5 and v = 2.72 x 10" 13 Hz s" 1 . 
The frequencies of maximum Y§ are plotted with this 
linear model subtracted if Y% > 35, with symbol areas 
proportional to Yq — 35. The abscissa is the observa- 
tion times folded with a 242 day period, with epoch 
MJD 48362 (the first day of useful BATSE data). A 
242 day period was used rather than the 241 ± 1 day 
orbital period estimate of Zhang et al. (1996) be- 
cause it resulted in less scatter. The horizontal bars 
give the four day width of the search intervals. The 



rapid frequency decrease at the beginning of the out- 
burst is undoubtedly due to the Doppler signature of 
an eccentric binary orbit. The slower increase in the 
observed frequency that follows could be due to both 
orbital doppler shifts and intrinsic accretion-induced 
spin-up of the pulsar. 

2.2. Pulse Frequency Modeling 

To obtain a preliminary estimate of the system pa- 
rameters and any intrinsic spin-up occurring during 
the outbursts, the frequencies in Fig. 2 were fit with a 
model which included an independent spin frequency 
and spin frequency rate for each outburst, and a full 
set of binary orbital parameters. A new search of the 
four day intervals was then conducted, and new fre- 
quency estimates obtained. In this search the segment 
phase offsets in equation 4 were replaced with 

A0 m = ^model(tm) + Aj/(t m - t) - <f>o(t m ) (9) 

where <ft mo de\(t) encorporates the spin- frequencies, fre- 
quency rates, and orbital parameters determined by 
the fit to the initial frequency estimates. This itera- 
tive step was needed to improve the sensitivity in in- 
tervals where rapid frequency changes were inferred, 
such as near the beginning of the outbursts. The 
search was restricted to frequency offsets less than 
5.8 /xHz from the new folding model. The model pa- 
rameters were then re-estimated using the barycentric 
frequencies determined for detections with Y" 6 > 35. 

The first column of Table 1 shows the estimated bi- 
nary orbital parameters, while the spin- up rates dur- 
ing the outbursts, and those inferred between out- 
bursts (using a 30 day outburst width) are given in 
Fig. 3. The highly eccentric orbit has a mass function 
of 7 ± 1M Q , suggesting a 9M Q companion or larger. 
The outbursts occur close in phase to the estimated 
periastron passage. The fit suggests that the long- 
term spin-up trend is the result of accretion torques 
that occur during the outbursts. 

The quality of this fit is poor, with \t = 34.3. The 
fit residuals for most of the outbursts show significant 
structure, indicating the simple constant frequency 
rate model used is inadequate. The errors given in 
Table 1 have been increased by the factor ^/xj to 
account for this poor fit, however they are undoubt- 
edly still underestimated. Adequate modeling will re- 
quire parameters describing the time history of the 
torque during the outburst. These parameters will be 
strongly coupled to the orbital elements due to the 



5 



Table 1 

Orbital Parameter Estimates 





Frequency Fit 


Phase Fit 


Period 


242.18±0.02 days 


242.180±0.012 days 


Periastron epoch 


JD 2449617.3±0.2 


JD 2449616.98±0.18 


Semi-major axis 


725±40 lt-s 


689±38 lt-s 


Eccentricity 


0.87±0.01 


0.8792±0.0054 


Argument of periastron 


264° ± 7° 


252.2° ±9.4° 




926/27 


55.1/52 



limited orbital phase coverage. Such modeling will be 
studied in the pulse timing analysis that follows. 

While statistically a poor fit, this model is suffi- 
ciently accurate for coherent epoch folding over in- 
tervals of a few days. The R.M.S. deviation of the 
measured frequencies from the model is 0.21 fiHz. 

2.3. Pulse Profile Measurements 

Fig. 4 shows the mean pulse profile from the most 
significant four day interval (MJD 49624.0-49628.0). 
The profile is derived from estimates of its first 6 
Fourier amplitudes, using the frequency model dis- 
cussed in section 2.2, and is therefore a smooth curve. 
Errors are shown at phases which are approximately 
statistically independent. The pulse amplitude is «20 
mCrab R.M.S. (20-50 keV). This profile shows a nar- 
row notch, which is similar to the GS 1843-02 profile 
observed by Ginga in the 9-38 keV range (Koyama et 
al. 1990a). 

In preparation for a pulse timing analysis, one 
mean pulse profile was estimated every two days dur- 
ing the time range of each outburst for which pul- 
sations were detected. Two profiles were made for 
each four day frequency search interval to provide 
higher resolution during the periastron passage. We 
found that profiles from one day integrations were of- 
ten not of sufficient significance to reliably estimate 
pulse phases. The first 6 Fourier amplitudes of each 
mean pulse profile were estimated from the 300s pro- 
files using the phase model obtained from the fit to 
pulse frequencies discussed above. The Poisson statis- 
tical errors of each Fourier amplitude [ik was modified 

by the factor (J h (n h )/[2M - 2]) / 2 , were M is the 
number of 300s profiles, to account for non-Poisson 



noise. Several profiles from the edges of outbursts 
were dropped from further analysis because of low 
significance. 

2.4. Pulse Phase Measurements 

A phase offset A<p from the phase model was esti- 
mated for each two-day mean pulse profile by fitting 
the Fourier amplitudes fih of the profile to the Fourier 
amplitudes of a scaled and shifted template profile: 

h=l a »h 

Here A is the estimated pulse amplitude, Th is a 
Fourier amplitude of the template profile, and a^ h 
is the error on the real or imaginary component of 

The template profile was obtained in three steps. 
Initially a cosine template was used to obtain phase 
offsets in the bright interval MJD 49620.0-49630.0. 
The profiles in this interval were then aligned using 
these offsets and a new template constructed from 
the weighted average of the aligned profiles. Using 
this template, phase offsets and amplitudes were than 
estimated for all the profiles. The 46 profiles with 
pulse amplitude above 5cr significance where than 
aligned using these estimated offsets, and the final 
template constructed from the weighted average of 
these aligned profiles. 

For each two-day interval j the mean observation 
time, reduced to the solar system barycenter, tj, was 
computed. The total pulse phase at that time was 
then computed as <j)j — 4> mo dci (tj ) + A<f)j , where Acfij 
is the estimated phase offset in interval j. 
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2.5. Pulse Phase Modeling 

The pulse phases were fitted with a model that as- 
sumed a Gaussian profile for the spin-up rate during 
each outbursts. In outburst k the spin-up rate is de- 
scribed by three adjustable parameters; the peak spin- 
up rate Ok, time of peak i^ cak , and outburst width 
Wfe. During outburst k the fitted phase model had 
the form 



, • peak rt pt 



del(*) - <Pk + "kit™ ~ T k ) 
h 



(11) 



ik exp 



dtdt 



Here the phase constant tpk and the frequency con- 
stant Vk are adjustable parameters, while the epoch 
Tfc is fixed near the center of the outburst. The emis- 
sion time i om is related to the barycenter corrected 
observation time t by 

t = t cm + a x sin i (sin u[cosE - e] (12) 
+ \/l — e 2 cos oj sin E)/c 



where 



E — e sin E 



2tt 



P. 



orbit 



Tpcriastron 



(13) 



with the orbit's projected semi-major axis a x sini, ar- 
gument of periastron uj, eccentricity e, orbital period 
P or bit, and periastron epoch Tp Crias tron, all being ad- 
justed in the fit. The fitting was performed by mini- 
mizing 



N 



X 



[<f>k ~ 



^moclcl 



(14) 



using the Levenberg-Marquardt method (Press et al. 
1992). 

Initially only a single common value was estimated 
for the widths Wk, and the times of peak tu were con- 
strained to be separated by the orbital period. This 
fit had a x 2 of 135.4 with 70 degrees of freedom. Al- 
lowing the Tfc to be independently adjusted reduced 
the x 2 1° 74.7 with 61 degrees of freedom. With the 
widths Wk also independently adjusted, x 2 dropped 
to 55.1 with 52 degrees of freedom. 

The orbital parameters obtained from this fit are 
shown in the second column of Table 1. These are 
consistent with those obtained by fitting observed 
frequencies. The period, periastron epoch, and ec- 
centricity are all well constrained by the pulse phase 
measurements. The argument of periastron uu and the 



semi-major axis a^sinz are not as well constrained. 
Fig. 5 shows their joint parameter confidence region. 
The broad minimum of x 2 for these two parameters is 
due to the difficulty in separating orbital and intrin- 
sic torque signatures with data restricted to a small 
interval of orbital phase. 

The estimated spin-up rate model is shown in Fig. 
6 which also shows the mean spin-up rate inferred 
between outbursts. Torques between outbursts con- 
tributing an estimated 5 ± 23% of the long term fre- 
quency change. The error on this estimate is coupled 
to those of u> and a x sini, with the estimated fraction 
varying from -0.33 to 0.26 for fits with fixed values 
of uj and a x sin i within the 50% confidence contour 
shown in Fig. 5. 

Fig. 7 shows the estimated delay of the spin-up 
rate peak from periastron passage. The mean delay is 
9.4±0.7 days, which is shown in the figure. The error 
on the mean delay is coupled to those of uj and a x sini, 
with the estimated mean delay varying from 8.1 to 
10.0 days for fits with fixed values of uj and a x sini 
within the 50% confidence contour shown in Fig. 5. 
The errors on the individual delays are correlated due 
to this coupling between orbital and torque parameter 
errors. In order to show the significant differences 
between delays, the error shown for each delay is that 
of the difference of the delay from the mean delay. 

Fig. 8 shows the estimated width parameters Wk 
of the spin-up rate profiles. This has a mean value 
6.2 ± 1.0 days, which is shown in the figure. The 
error on the mean width is coupled to those of uj and 
a x sini, with the estimated mean varying from 5.2 to 
7.6 days for fits with fixed values of ui and a x sini 
within the 50% confidence contour shown in Fig. 5. 
The errors on the individual widths are correlated due 
to this coupling between orbital and torque parameter 
errors. In order to show the significant differences 
between widths, the error shown for each width is that 
of the difference of the width from the mean width. 

2.6. Outburst Profile 

Fig. 9 shows the lightcurve of the pulsed flux for 
the 1996 September outburst. The other outbursts 
are similar. A mean pulse profile was determined from 
two day intervals, using the detailed phase model es- 
timated in the previous section to align the 300s seg- 
ment profiles. Non-Poisson noise was accounted for 
in the errors as previously described. The pulsed flux 
for each interval was calculated from the correlation 
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of the mean pulse profile in that interval with the 
pulse template (with no adjustment in phase). The 
template was normalized to have unit variance, so 
that the correlation gives the R.M.S. pulsed flux if 
the pulse profile has the same shape as the template. 
Also shown in the plot is the model torque profile for 
the outburst, and the time of periastron passage. 

To examine the outburst lightcurve with higher sig- 
nifigance an average outburst profile was constructed. 
This is shown in Fig. 10. The pulsed flux was deter- 
mined as above for P or bit/242 « 1 day intervals and 
then epoch folded at the orbital period P or bit- Also 
shown in the figure is the average model torque pro- 
file. Overall the pulsed flux and the torque profile are 
well aligned and of approximately the same width. 
Since the widths and peak times of the torque profile 
for each outburst were estimated in the phase model- 
ing, this demonstrates a correlation between the spin- 
up rate and the pulsed flux. 

There is significant narrow feature in the average 
outburst lightcurve near periastron passage. This has 
a width of two days or less, and centered approxi- 
mately a half day before periastron. Inspection of the 
10 individual outburst lightcurves at 1 day resolution 
shows that this feature is present in the majority of 
them, and is not caused by any individual flux en- 
hancement. 

In figure 11 we show the average total flux pro- 
file obtained by epoch folding, at the orbital period, 
fluxes obtained by Earth occultation analysis (Har- 
mon et al. 1992) of the BATSE CONT data. This 
figure includes only the outbursts analysed by Zhang 
et al. (1996) using data from MJD 48363-49741. The 
occultation fluxes were obtained from fits that used a 
power-law spectral model with a fixed photon num- 
ber index of -2.5. The onset time and width of the 
outburst roughly agrees with that in figure 10, but a 
detailed comparison is precluded by the higher rela- 
tive noise level in the average total flux profile. 

3. Discussion 

These observations show a series of 10 outbursts of 
the pulsar GS 1843-02 occurring every 242 days. This 
outbursting behavior, and the wide, highly eccentric 
orbit, which has a mass function requiring a compan- 
ion mass of M opt > 7M Q , leave little doubt that the 
source is a Be star/pulsar system. The consistency 
of the outburst intervals with those of the periodic 
transient GRO J1849-03 (Zhang et al. 1996) clearly 



shows these to be the same source. BeppoSAX Wide 
Field Camera observations of GRO 1849-03 during a 
predicted outburst in September 1996 have recently 
identified this source with 2S 1845-024 (Soffitta et al. 
1998). 2S 1845-024 was observed by SAS 3 in July 
and September 1975, and localized to 35" (Doxsey 
et al. 1977). As shown in Soffitta et al. (1998) the 
recent outbursts were also seen in the RXTE ASM 
lightcurve of 2S 1845-024. 

Other observations of outbursts and additional 
identifications might possibly be found by examin- 
ing sources cataloged in the region of 2S 1845-024 by 
early missions. This includes A 1845-02 and H 1845- 
024. A 1845-02, which was discovered by Ariel V in 
1974 November (Villa et al. 1976), and was identified 
with 2S 1845-024 by Doxsey et al. (1977). H 1845-024 
is listed in the HEAO 1 catalog (Wood et al. 1984) 
with the assumed location of 2S 1845-024. In Table 
2 we list observations of 2S 1845-024 and its possi- 
ble counterparts along with the calculated orbit phase 
during the observations. The detections by Ariel V, 
SAS 3, Ginga, BeppoSAX, and RXTE are all consis- 
tent with the range of orbital phases in which BATSE 
detects pulsations, which is from 4 days before to 24 
days after periastron passage. The HEAO observa- 
tion however occurs near apastron. This detection is 
possibly due to source confusion with another tran- 
sient. Several other transient pulsars that have been 
found in this region (Koyama et al. 1990b), since the 
HEAO 1 catalog analysis. 

BATSE observations of accreting pulsars (Bild- 
sten et al. 1997) have shown that extended series 
of low luminosity (10 36 ~ 37 erg s _1 ) outbursts occur- 
ring near periastron are common in the transient Be 
star/pulsar systems. The outbursts of 2S 1845-024 
observed by BATSE fit into this "normal" outburst- 
ing pattern. The peak luminosity determined for 
the 1996 September outburst from combined RXTE 
ASM (2-12 keV) and BATSE occultation data (20- 
100 keV) was 6 x 10 36 erg s -1 , assuming a distance 
of 10 kpc (Soffitta et al. 1998). As shown in Ta- 
ble 2 the peak fluxes of the outbursts observed by 
BATSE are all fairly similar. The outbursts have 
proceeded like clockwork since at least the launch 
of the Compton Observatory. Fig. 12 shows the 
the Ginga and BATSE frequency measurements, cor- 
rected for the orbit in the second column of Table 1. 
Significant spin-up occurred between the Ginga and 
and the BATSE measurements, but that spin-up was 
smaller than would be predicted from an extrapola- 
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tion of the BATSE measurements. If, as suggested by 
our estimates, the long-term spin-up is due largely to 
torques occurring within outbursts, then unobserved 
outbursts must have occurred between the Ginga and 
BATSE observations. The lower spin-up rate implies 
that these outbursts were on average smaller, or that 
one or more periastron passages occurred without an 
outburst. 

The BATSE observations of 2S 1845-024 provide 
the best evidence to date of spin-up occuring dur- 
ing normal outbursts of a transient pulsar. Spin-up 
during the less common "giant" outbursts of tran- 
sient pulsars has been reported many times. During 
giant outbursts, which peak near the Eddington lu- 
minosity, an accretion disk is known to be present 
(Finger, Wilson & Harmon 1996a). Only recently, 
however, has evidence for spin-up during normal out- 
bursts emerged with BATSE observations of 2S 1417- 
624 (Finger, Wilson & Chakrabarty), GS 0834-430 
(Wilson, C. A., et al. 1997), and EXO 2030+375 
(Stollberg et al. 1999), and now the observations 
of 2S 1845-024 present here. The most common ex- 
planation for normal outbursting activity in transient 
Be star/pulsar systems is the direct accretion of the 
dense equatorial wind (or circumstellar disk) of the Be 
star (Waters et al. 1989) as the neutron star passes 
through periastron. Yet wind accretion is thought to 
be very inefficient at transfering angular momentum 
(Ruffert 1997), leading us to expect any spin- up dur- 
ing normal outbursts to be fairly small. 

The spin-up rate we see during outbursts of 2S 
1845-024, and its correlation with the pulsed flux 
(Fig. 10) suggests an accretion disk is present. A 
disk will form if the specific angular momentum I 
in a wind exceeds the Keplerian specific angular mo- 
mentum l m — (GMr m ) 1 2 at the magnetospheric ra- 
dius r m . After an accretion disk is formed, the spe- 
cific angular momentum of accreting material is main- 
tained near l m . For accretion to occur, the magneto- 
spheric radius must be within the co-rotation radius 

1/ _ 2/ 

r co = (GM) 1 3 (2itv) 1 3 . Therefore the specific an- 
gular momentum needed for disk formation, l m , is 
bound by 

l m < lco = (GMr co ) V2 = 8 x lO^cmV 1 . (15) 

From the luminosity and spin-up rate at the peak of 
the 1996 September outburst (which is typical) we 
infer a specific angular momentum of the accreting 



material of 

I = ItxIvM- 1 = 8 x 10 17 di" 2 cmV 1 (16) 

assuming M = L(GM/R)-\ Here v = 4xl0~ 12 Hzs" 1 
is the spin-up rate, M is the mass accretion rate, 
L = 6 x 10 36 erg s _1 d\ is the luminosity, dm is 
the distance in units of 10 kpc. We have assumed 
I = 10 45 g cm 2 for the neutron star moment of iner- 
tia, M — 1.4M Q for the mass, and R = 10 6 cm for 
the radius. From comparison of equations 15 and 16 
we conclude that disk accretion is occuring during 
the outbursts of 2S 1845-024. By moving the source 
farther away we could avoid this conclusion. The dis- 
tance of 10 kpc estimated by Koyama et al. (1990a) 
was, however, based on the minimum column density 
determined for a set of Ginga spectra. Since part of 
this column density may be intrinsic to the source, a 
closer distance is more likely than a farther one. 

The progenitor system of a Be/neutron star binary 
is likely a compact binary system of B stars (van 
den Heuvel 1994). The more massive star evolves 
more rapidly and transfers mass to its companion via 
Roche-lobe overflow, resulting in a helium star in a 
circular orbit about its mass inhanced companion. 
The helium star then undergoes a supernova explo- 
sion producing a neutron star in a wide eccentric or- 
bit, if the system remains bound. Because the kick 
velocities produced by the supernova are thought to 
be larger than or comparable to typical orbital veloci- 
ties of the pre-supernova systems, many such systems 
are expected to be disrupted. 

Fig. 13 shows the kick velocity required for vari- 
ous possible progenitor systems to have produced the 
2S 1845-024 system (see e.g. Kalogera et al. (1996) 
for the relevant equations). Before the supernova the 
orbital angular momentum and the companion's spin 
angular momentum are aligned, while afterward they 
are separated by the angle 0. The calculation as- 
sumes a neutron star mass of 1.4M Q , a pre-supernova 
helium star mass of 3.0M Q , and a companion star 
mass of 11.3M Q (i = 60°). The current eccentricity, 
period, and semi-major axis are assumed unchanged 
since the supernova. These elements are unaffected by 
the tidal interactions and spin-orbital coupling which 
lead to periastron advance and precession of the or- 
bital plane and the Be star spin axis (Lai, Bildsten, & 
Kaspi 1995). For these parameters a minimum kick 
of 41 km s _1 is required, which occurs in a 595 day 
period system with kick directed along the direction 
of motion. On the other hand, a kick of 550 km s _1 
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Table 2 

Orbital Phase of Observations 



Date 


Phase 


Peak Flux 




(days past periastron) 


(mCrab) 


Ariel V a 




(2-6 keV) 


1974 Nov. 12 - 25 


11 - 26 


12 


1975 July 21 - Aug. 16 


21 - 48 


2 


SAS 3 b 




(2-11 keV) 


1975 July 11-16 


11 - 17 


9 


HEAO l c 




(2-10 keV) 


1977 Oct. 2-10 


98 - 106 


7 


Ginga LAC d 




(2-20 keV) 


1987 Oct. 6 


-120 


< 0.25 


1987 Oct. 10 


-116 


< 0.25 


1988 Sept. 12 


-20 


< 0.25 


1988 Sept. 29 


-3 


8 


1988 Oct. 3 


1 


9 


BATSE 




(pulsed, 20-50 keV) 


1991 May 28 - June 16 


-2 - 18 


18 ±5 


1992 Jan. 25 - Feb. 13 


-2 - 18 


12 ±2 


1992 Sept. 21 - Oct. 12 


-4 - 18 


14 ±2 


1993 May 23 - June 11 


-2 - 18 


15 ±2 


1994 Jan. 18 - Feb. 14 


-4 - 24 


15 ±2 


1994 Sept. 21 - Oct. 12 


-1 - 22 


14 ±2 


1995 May 19 - June 7 


-3 - 17 


20 ±4 


1996 Jan. 16 - Feb. 8 


-3 - 21 


13 ±2 


1996 Sept. 16 - Oct. 5 


-1 - 19 


14 ±3 


1997 May 14 - June 8 


-3 - 23 


15 ±2 


BeppoSAX Wide Field Camera 6 


(5-21 keV) 


1996 Sept. 17-18 


- 1 


9 


Rossi X-ray Timing Explorer All Sky Monitor 


(2-12 keV) 


1996 Sept. 22 - Oct. 3 


5-17 


21 ±2 


1997 May 20 - May 27 


3-11 


22 ±2 



a A 1845-02 was detected in scans folded over the given time ranges (Villa 
ct al. 1976). Fluxes were calculated from the HE count rates. Cross scans 
are mentioned, but accurate dates are not given. 

b From Doxsey et al. (1977). 

c From Wood et al. (1984). The observation interval was calculated from 
the ecliptic longitude of 2S 1845-024. 

d From Koyama et al. (1990a). The fluxes were calculated from the 2-20 
keV count rates. 

°From Soffitta et al. (1998). 

f From 4 day averages using 1 Crab = 75 cs _1 . 
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directed opposite to the direction of motion in a 9.7 
day period system could have produced the system. 
In this case the Be star spin and orbital angular mo- 
mentum are in opposite directions. The actual origin 
of the current system probably lies between these ex- 
tremes. In particular we wish to emphasis that if 
typical kick velocities are ~ 250 km s _1 (Hansen & 
Phinncy 1997), then there is no reason to suppose 
that the orbital plane and the Be star equator are 
closely aligned, as is normally assumed in the direct 
wind accretion outburst model (Waters et al. 1989). 

We propose that the pulsar's orbit is inclined to its 
companion's equator, and that the narrow feature in 
the outburst profile (Fig. 10) near the onset of the 
outburst is due to direct wind accretion during the 
pulsar's passage through the companion's circumstel- 
lar envelope. Fig. 14 shows the pulsars separation 
from its companion, and its true anomaly, for times 
close to periastron. The geometry is changing rapidly, 
with the separation changing by a factor of two in less 
than four days, and the true anomaly changing by 43° 
per day at periastron. If the angle between perias- 
tron and the node of the orbit with the companion's 
equator is less than 40°, then passage through the 
equatorial plane will occur within a day of periastron 
passage. Be star circumstellar envelopes have half- 
angles of < 10° (see Hanuschik 1996, Quirrenbach et 
al. 1997, Wood et al. 1997). The passage through the 
circumstellar envelope would therefore take less than 
a day for orbits inclined from the companion's equa- 
tor by more than 15°. Direct wind accretion during 
the passage could result in the onset feature in the 
profile. At this time the accretion disk could also 
be forming. Alternatively a disk could already exist 
(Bildstcn et al. 1997), with enhanced accretion initi- 
ated near periastron due to tidal torques and newly 
added material. 

This neutron star's wide and eccentric orbit about 
a massive star is similar to two radio pulsars PSR 
B1259-63 (Johnston et al. 1992) and PSR J0045- 
7319 (Kaspi et al. 1994). Both of these pulsars are 
in highly eccentric orbits (e = 0.8698, 0.8080) about 
8 — 10M Q main sequence stars with periods of 1236 
and 51.2 days respectively. Bildsten et al. (1997) 
had earlier noted the lack of accreting pulsars with 
similarly large eccentricities, and conjectured it was 
a selection effect due to the difficulties of repeatedly 
observing long period systems. 2S 1845-024 begins to 
fill in this gap. As we have suggested for 2S 1845- 
024, the orbital plane of PSR J0045-7319 is believed 



to be inclinated from its companion's equator (Lai, 
Bildsten, & Kaspi 1995, Kaspi et al. 1996). However, 
given the its low spin frequency, it is unlikely that 2S 
1845-024 would be a radio pulsar now or in the future. 

M.H.F. acknowledges support from NASA grant 
NAG 5-4238. L.B. and T.A.P. acknowledge support 
from NASA grant NAGW-4517. We thank the anony- 
mous referee for helpful comments, which have im- 
proved the text. 
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Fig. 3. — The intrinsic spin-up rates estimated by fit- 
ting the pulse frequencies. The model assumes the 
spin-up rate is constant for 30 day intervals contain- 
ing each outburst. The rate between outbursts is cal- 
culated from the model frequency difference between 
the edges of the adjacent outburst intervals. 



This 2-column preprint was prepared with the AAS IATJrjX 
macros v4.0. 
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Fig. 4.— The pulse profile of GS 1843-02 obtained 
from the BATSE 20-50 keV rates. The profile is lim- 
ited to six harmonics. Error bars are give at approx- 
imately statistically independent points. 



Fig. 5. — Joint confidence regions for the argument 
of periastron to and the projected semi-major axis 
a x sinz as estimated by fitting pulse phases. 
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Fig. 6. — The intrinsic spin-up rates estimated by 
fitting the pulse phase. The model profile is Gaus- 
sian with adjustable width and peak time for each 
outburst. The solid circles give the peak spin-up 
rates. The average spin-up rates between outbursts 
are given by the solid diamonds. These are calcu- 
lated from the estimated model for 182 day intervals 
centered between outbursts. 
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Fig. 7. — The times of peak spin-up rate ^P eak relative 
to periastron passage estimated from the pulse phase 
fit. The errors shown are for the difference from the 
mean delay of 9.4±0.7 days, which is shown by the 
dashed line. 
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Fig. 8. — The estimated widths Wk of the spin-up 
rate profiles. The errors shown are for the difference 
from the mean of 6.2±1.0 days, which is shown by the 
dashed line. 
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Fig. 9.— The 20-50 keV pulsed flux (R.M.S. nor- 
malized) during the 1996 September outburst. Also 
shown is the model torque profile for this outburst 
(smooth curve). The dotted line is the time of peri- 
astron passage. 
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Fig. 10.— The 20-50 keV pulsed flux (R.M.S. normal- 
ized) epoch folded at the orbital period. Also shown 
is the average model torque profile. 
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Fig. 11— The 20-50 keV total flux (from BATSE 
Earth occultation measurements) epoch folded at the 
orbital period. 



10 





i 




i 

j 

l 

i 

t 

1 

1 




1 

f 

» 


s 





7000 



8000 9000 10000 

Time (JD-2440000.5) 



1 1000 



Fig. 12. — The long-term history of the spin fre- 
quency, corrected using the orbit in Table 2. The ini- 
tial points are from the Ginga pointed observations, 
the remainder from BATSE observations. 
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Fig. 13. — Kick velocity (km s _1 ) required to produce 
the present orbit for prc-supernova systems with a 
companion mass of 11.3Mq, and a neutron star mass 
of 1.4M Q . 9 is the angle between the orbital angular 
momentum and the Be star spin angular momentum 
after the supernova explosion. 
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Fig. 14. — The neutron star's separation from its com- 
panion and the true anomaly, for times near perias- 
tron passage. The separation is calculated assuming 
a neutron star mass of 1.4M©, and an inclination of 
60°. 



